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Abstract 

It is shown that the technique recently suggested by Lowell Brown for summing the 
tree graphs at threshold can be extended to calculate the loop effects. Explicit result 
is derived for the sum of one-loop graphs for the amplitude of threshold production 
of n on-mass-shell particles by one virtual in the unbroken A0 4 theory. It is also 
found that the tree-level amplitude of production of n particles by two incoming 
on-mass-shell particles vanishes at the threshold for n > 4. 



The problem of calculating amplitudes of processes with many weakly interact- 
ing particles has recently attracted a considerable interest, initially triggered by 
the observationE S that such processes in particular are associated with a pos- 
sible baryon and lepton number violation in high-energy electroweak interactions. 
Cornwall@ and Goldberg^ have pointed out that in perturbative amplitudes with 
many external particles the weak coupling may get compensated by large number of 
diagrams. This is a manifestation of the old-standing problem of the factorial growth 
of the coefficients in the perturbation theoryB . Since the perturbative expansion for 
multi-particle amplitudes starts from a high order in the coupling constant, for suffi- 
ciently large number n of particles the factorial growth of the coefficients in the series 
invalidates the perturbative calculation of such amplitudes. Given lack of a better 
approach it seems useful to quantify and study the problem within the perturbation 
theory itself. A simple model example in which the problem arises with full strength 
is the amplitude A n where a virtual particle of a real scalar field produces a large 
number n of on-mass-shell 0-particles in the A0 4 theory. It has been recently found 
that the sum of all tree graphs for this amplitude in the threshold limit, i.e. when all 
the produced particles are at rest, can be calculated exactly for arbitrary n both in 
the case of unbroken symmetry^ and in the case of theory with spontaneous breaking 
of the symmetry under the reflection — > — S. 

Originally the calculationS was done by directly solving a recursion relation for 
the tree graphs. Argyres, Kleiss and PapadopoulosB applied a regular method of 
solving the recursion relations based on generating function for which the recursion 
relation for the amplitudes A n is equivalent to second order non-linear differential 
equation. Most recently Brown0 has shown that the generating function is nothing 
else than the classical field <fio(t) generated by an external source p = poe imt which field 
is a complex solution of the Euler-Lagrange classical equation satisfying the condition 
that it has only the positive frequency part. The equations and their solutions in both 
approaches are related by a simple change of variable. Thus Brown has reproduced 
the previous results^ in a simple and elegant way. 

The purpose of this paper it to show that Brown's technique can be extended to 
calculate the loop contributions to the amplitudes A n as well as the amplitudes of 
more complicated processes e.g. of the scattering 2 — > n. For definiteness the case of 
unbroken A0 4 theory with the Lagrangian 
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will be considered and we will concentrate on the calculation of the one-loop correction 
to the threshold amplitudes A n . The result of this calculation can be written as 



(2fc + l|0(O)|O) = (2fc + l)! 
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where A and m are the renormalized coupling constant and the mass, the appropriate 
renormalization condition, which specifies the finite terms, will be described below. 
Equation (fj) gives the exact sum of all tree and one-loop graphs at the threshold 
of production of n = 2k + 1 ^-particles for arbitrary k []. The formula (0) gives 
the relative magnitude of the first loop correction growing as n 2 at large n. This 



behavior explicitly demonstrates invalidity of the previous arguments^ 1 III) of the 
present author that terms, containing n 2 X should be absent in the loop effects. The 
reason for those arguments to be faulty is related to the singularities of the underlying 
classical field in the plane of complex time, where the quantum fluctuations are more 
singular than the classical background. However, the presence of the n 2 A parameter 
in the loop effects does not necessarily imply that the quantum effects completely 
eliminate the growth of the amplitudes and a further study is needed. The relatively 
easy calculability of the threshold amplitudes at the tree and the one-loop level and 
the remarkably simple form of the result (fj) gives us a hope that these amplitudes 
can be studied well beyond the first terms in the perturbative expansion. 

That the amplitudes in the A0 4 theory at the multi-particle thresholds may have 
special properties is also hinted at by another fact, which follows as a by-product from 
the calculation in this paper. Namely, if one considers the sum of all tree graphs for 
the amplitude of the process where two incoming on-shell particles produce n on-shell 
particles exactly at the threshold, it turns out that this sum is non-zero only for n = 2 
and n = 4 and is vanishing for all n > 4 (the number of final particles in this case is 
necessarily even). It is in fact due to this behavior that the one-loop term in eq.(0) 
contains only the 4-particle threshold factor. 

The technique suggested by BrownS is based on the standard reduction formula 

lr The number n of final particles produced by one virtual is necessarily odd due to the unbroken 
reflection symmetry. 
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representation of the amplitude through the response of the system to an external 
source p(x), which enters the term pip added to the Lagrangian 



(n|0(x)|O) 



n . 

J] lim / d A x a e ipaXa (m 2 -p 2 a ) 

2=1 
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(O out |0(a;)|O m ) p | p= o , (3) 



the tree-level amplitude being generated by the response in the classical approxima- 
tion, i.e. by the classical solution <po(x) of the field equations in the presence of the 
source. 

For all the spatial momenta of the final particles equal to zero it is sufficient to 
consider the response to a spatially uniform time-dependent source pit) = p (uj)e lult 
and take the on- mass-shell limit in eq.@ by tending uj to m. The spatial integrals 
in eq.@ then give the usual factors with the normalization spatial volume, which as 
usual is set to one, while the time dependence on one common frequency u implies 
that the propagator factors and the functional derivatives enter in the combination 



(m 2 -p 2 ) — — - -> (m 2 - uj 2 )— — - = — — - , (4) 
1 Fa! 5p{x a ) y J 5p(t) 5z(t) V ; 



where 



Po (u)e^ 

Z(t) = — ; 5 

m z — xe — uj z 

coincides with the response of the field to the external source in the limit of absence 
of the interaction, i.e. of A = 0. For a finite amplitude p of the source the response 
z(t) is singular in the limit uo — > m. The crucial observation of BrownS is that, since 
according to eq. (^) we need the dependence of the response of the interacting field <f> 
only in terms of z(t), one can take the limit Pq{uj) — > simultaneously with to — > m 
in such a way that z(t) is finite, 

z(t)^z e mt . (6) 

Furthermore, to find the classical solution <po(x) in this limit one does not have 
to go through this limiting procedure, but rather consider directly the on-shell limit 
with vanishing source. The field equation with zero source is of course given by 

<9 2 + m 2 + A0 3 = . (7) 
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For the purpose of calculating the matrix element in eq.(|3|) at the threshold one 
looks for a solution of this equation which depends only on time and contains only 
the positive frequency part with all harmonics being multiples of e tmt . The solution 
satisfying these conditions reads as®: 

A, U\ Z ® 

<pa(t) 



1 - (X/8m 2 )z(t) 2 

According to equations (£|) and ([3D the n-th derivative of this solution with respect 
to z gives the matrix element (n|0(O)|O) at the threshold in the tree approximation: 

(2* + l|0(O)|O>o=^J 0oU = (2A; + l)!^j , (9) 

which reproduces the previously known result^. The fact that the matrix element 
is non-zero only for odd n obviously follows from that the expansion of 0o in e( l-@ 
contains only odd powers of z. 

It can be noticed that the solution (H) is in fact not uniquely determined by the 
above mentioned conditions. Namely, z(t) can be rescaled by an arbitrary constant 
C . This constant corresponds to the choice of normalization of the field, so that the 
value C — 1 is fixed by the usual normalization condition (1|0(O)|O) = 1, as can be 
seen from the linear term in the expansion of 4>o in powers of z. 

Another important point concerning the solution (Q) is related to the fact that 
this solution is essentially complex for real time t. This is imposed by the that in 
calculating production of particles by the virtual field, rather than both production 
and absorption, one necessarily has to consider only the positive frequency part of 
the field, which is essentially complex. 

The quantum loop corrections to the amplitudes (n|0|O) are obtained by substi- 
tuting instead of the classical field the mean value of the full field 

(j>{x) = (j>o{x) + <j>g{x) , (10) 

where <f> q (x) is the quantum part of the field. Expanding the field equation (|7|) near the 
classical solution O an d retaining only the first non-vanishing quantum correction, 
one finds that the mean field <p(x) to the first quantum order satisfies the equation 

d 2 <f){x) +m 2 0(x) + A0(x) 3 + 3A0 o (x)(0 9 (x)0 g (x)) = , (11) 
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where ((j) q (x)(f) q (x)) is the limit of the Green function in the classical background field 

00 



G(x 1 ,x 2 ) = (T(<j> q ( Xl )<l> 9 (x 2 ))) (12) 

when its arguments are at the same point x. 

Therefore the steps needed to calculate the first loop correction to the amplitudes 
A n are the following: 

i. Calculate the Green function (|12|) as the inverse of the operator of the second 
variation of the action, 

d 2 + m 2 + 3\<po(x) 2 ; (13) 

ii. Find its limit in coinciding points, which enters eq. flTl"]) ; 

Hi Expand the solution of thus found equation in powers of z(t), which gives the 
amplitudes at the threshold in the same way as in eq.([|). 

This program however is obscured at the very first step by the fact that with es- 
sentially complex 0o ( e Q-@) the operator (|T3"| ) is essentially non-Hermitean. However 
one can render this operator real and thus the problem more tractable by analytical 
continuation in time t, which amounts to rotation and shift in the complex plane. 
Namely, the substitution which achieves the goal reads as 



A z{t) = ie mT (14) 



T = it + — m^-| - — (15) 



and the variable 

1 A^o in 
— In 

m 8m 2 2m 

is then used as the new time variable t. The necessity of the shift in addition to the 
usual rotation to the Euclidean time is caused by existence of a pole of <fio(t) on the 
negative imaginary axis, where the operator (|i~3l) is singular. The poles are repeated 
parallel to the real axis with the period n/m. The axis, corresponding to real r, 
on which the operator (|i~3l) is real runs parallel to the imaginary axis of t exactly 
in the middle between two poles, see Fig.l. It should be emphasized however that 
it is the pole structure of the field, which gives rise to the factorial growth of the 
multi-particle amplitudes the study of which may eventually be the central point in 
solving the problem of multi-boson processes. Here, for the purpose of the specific 



5 



calculation, we chose to stay away from the poles to avoid explicit singularity in the 
equations. 

To somewhat simplify the notation we set the mass m equal to one and restore it 



when needed and also introduce the notation u(r) = e T = — iJ\/8z(t). For real u(r) 
the classical field (§) is purely imaginary: 

= i/n^VissT (16) 

and the operator ( |i~3l) is real. In a mode with spatial momentum k the operator has 
the form 

d 2 2 6 
~d^ + LJ ~ (cosh rf ' (17) 

which is the familiar operator in one of exactly solvable potentials in Quantum Me- 
chanics (see e.g. Ref.|[L2|), and u is the energy of the mode: u 2 = k 2 + 1. 

The regular at r — > +oo solution of the homogeneous equation with the operator 
( |T7D has the form 

» / / % % 2-3u + u 2 -8m 2 + 2uo 2 u 2 + 2u 4 + 3uju 4 + uj 2 u 4 

fMr)) = — — 2 (18) 

(1 + u z ) 

and the solution regular at r — > — oo is given by 



, , - » {2 + Zlo + lo 2 -%u 2 + 2u 2 u 2 + 2u 4 -Zuu 4 + u 2 u 4 ) 

f2{u{r)) = fi{l/u{r)) = 2 . 

(1 + u ) 

(19) 

The Wronskian of these solutions is given by 

W = h{r)f 2 {r) - f[{r)f 2 {r) = 2co(co 2 - l)(co 2 - 4) . (20) 

The convention for the sign of the Green function used here is specified by the explicit 
expression for the Green function in partial wave with the spatial momentum k in 
terms of f h f 2 and W: G uj (t 1 ,t 2 ) = fi(n)f 2 (T 2 )/W for n > r 2 and G w (t x ,t 2 ) = 
h{r 2 )f2{ri)/W for r 2 > n. 

Naturally, having the explicit expression for the Green function one can also eval- 
uate amplitudes of more complicated processes, say, the tree level amplitude of the 
threshold production of n particles by two incoming on-mass-shell particles of high 
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energy. However equations ( |18D , fll~9|) and fl2"0|) show that there is in fact almost noth- 
ing to calculate for the latter amplitude: the Green function has poles only at uj 2 = 1 
and uj 2 = 4 (the zeros of the Wronskian (pOD). By the reduction formula this implies 
that the on-mass-shell amplitude is non-vanishing only at these values of the energy 
of each of the two incoming particles. The case uj = 1 corresponds to the trivial 
process 2 — > 2 at the threshold, while the case uj = 2 corresponds to the process 
2^4. (In the rest frame of the produced particles, which is used throughout this 
paper, u corresponds to the energy of each of the two incoming particles, so that the 
total energy is 2ou — 4.) The absence of other poles of the Green function at higher 
uj means that for n > 4 the sum of tree graphs for the on-shell process 2 — > n is 
vanishing. 

After this remark we proceed with calculating the loop correction in eq.(0). The 
partial wave Green function at coinciding points is given by 

g U (T) = G u {T,T) = f 1 {T)f i (T)/W (21) 

which yields the average value of the square of quantum fluctuations after integrating 
over k: 

d 3 k 



(Mr)Mr)} = 9(r) = J ft-OO^a = ^ ]j u {t)u>V<J=I duj . (22) 

The calculation of the latter integral involves problems related to the on-shell singu- 
larities and to the ultra violet divergence. The on-shell singularities correspond to 
the zeros of the Wronskian (^) at uj 2 = 1 and uj 2 = 4. The first of these corresponds 
to the translational zero mode of the classical solution <fio and in fact produces no 
effect in the integral in eq.(|22|) since the singularity at uj 2 = 1 is integrable. (This 
is why in the 4-dimensional theory one does not have to consider subtraction of the 
contribution of the zero mode from the Green function.) The pole at uj 2 = 4m 2 (the 
dependence on mass is restored) is dealt with using the Feynman's ie rule i.e. by 
shifting the pole to the negative half-plane m 2 — > m 2 — ie. The integral then develops 
imaginary part, which in the end corresponds to the dynamical imaginary part of the 
one-loop graphs, dictated by the unitarity. 

To separate the ultra violet divergent terms we expand g w {r) in powers of uo~ l 
and find that the two terms, which give the quadratic and the logarithmic divergence 
have the form 
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1 6m 2 

S "< r) " + (TT^)V +9 " (T) ' (23) 

where the regular part <j£(t) contains terms of the order u~ b and higher, so that 
its contribution to the integral in eq.([22|) is finite in the ultra violet. After this 
decomposition the result of the integration in eq. ( p2]) can be presented as 

9(r) = ~/i + + A) - tt^-^F , (24) 



where 



2vr 2 V 2 - y/3 
and Ji, I3 are the ultra violet divergent integrals: 



2 * 1 + 2V " 27T 2 ' 1 + M 



^^i±^_ !T ) ( 25 ) 



2tt 2 



u 1 ^ 



Vuj 2 - 1 dcj . (26) 



In equation (|24|) we have also combined with the logarithmically divergent integral 
a part of the finite contribution from integration of the 5^(t), which has the same 
functional dependence on u(r), hence the factor (I 3 + l/(27r 2 )). 

The divergent contributions can be regularized in a standard way, the most 
straightforward being the Pauli-Villars regularization. Upon substitution into equa- 
tion (|ll]) for the mean field with the quantum correction the quadratically divergent 
part proportional to I\ gives rise to a term linear in the classical field 0o while the log- 
arithmically divergent part proportional to -Z3 results in a correction to the term with 
A0Q- Therefore these terms can be dumped into the definition of the renormalized 
mass m and the coupling constant A according to 



m 2 = m 2 H 7i 

2 

9A 2 1 

A = A -- (/ ' + 2^»- (27) 

These definitions can be used to relate the quantities fh and A to the renormalized 
constants in any other renormalization scheme. One can readily see that the divergent 
parts are scheme-independent, while the relation between the finite parts depends on 
the specific definition of the regularization procedure. 
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The only non-trivial modification of the average field given by equation ((TTJ) is 
related to the finite part of the average value of the square of quantum fluctuations 
(eq.(^)), proportional to the constant factor F. If one seeks the solution of the 
equation (O) in the form = <pQ(t;rh,\) + <fii(t), where the renormalization of 
the constants is plugged into the functional dependence of the classical solution, the 
equation for the correction 0i(r) (i.e. on the r axis) reads as 

( d 2 24m 2 \ [~8 

(s? " 1 + (TTW ) 01 = \ \ JTTuPf ' (28) 

the condition on the appropriate solution to this equation being that its expansion in 
u starts with the fifth power, since only starting from final states with five particles the 
threshold amplitudes develop an imaginary part, which in this calculation originates 
in the imaginary part of F. The solution satisfying this condition is 

, , n 3A /8~ u 5 . , 

Mr) = ~ ! TV a (TTii 3 ) 1 ' (29) 

. Using equation (|T4| ) one can readily restore from here the response of the field in 
terms of z(t) with the first quantum correction included: 

6 m- *(*) U 3A F (V 8 ™ 2 ) 2 ^*) 4 ^ m) 

^ 0+l[ ) 1 - (A/8m 2 )^) 2 { 4 (1 - (A/Sm 2 )^) 2 ) 2 J 1 j 

and by expanding in series in z(t) finally arrive at the result in equation (0). 

The rotation (0) used here may invite the objection, that such rotation in the 
path integral is obstructed by the infinite chain of poles parallel to the real axis of 
t, which may give rise to extra contributions in the quantum effects. However it can 
be explicitly shown that this does not happen at least at the one-loop level. Namely, 
it is a straightforward (but rather cumbersome) exercise to verify that the recursion 
relations for the sum of graphs for the propagator of the field (j) with emission of n 
on-shell particles all being at rest are equivalent to the differential equation for the 
Green function of the operator fllTD and then that the recursion relations for the loop 
graphs are equivalent to the equation (|28|) on the r axis. Another simple (and in no 
way rigorous) check is to verify the formula (0) for few first n by direct computation 
of the graphs. This also turned to be helpful in checking the relative coefficients and 
signs in the equations of this paper. The remarkably simple form of the result (0) 
suggests that there may be a way to calculate further quantum effects. In particular 
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one can notice that the finite term, proportional to the factor F (eg. (|25D) has the 
form given by the simple scalar vacuum polarization at q 2 = 16m 2 . This of course is 
a consequence of the eigenmode of the operator (|T7|) at oj 2 = 4, or, equivalently, of 
the fact that the tree-level threshold amplitudes of the processes 2 — > n are equal to 
zero for n > 4. In terms of the graphs the cancellation of the contributions to the 
imaginary and the real parts of the thresholds at higher q 2 looks quite surprising. 

In the present calculation we have avoided approaching the poles of the classical 
solution <fto, where the quantum expansion in fact breaks down, since the quantum 
fluctuations are more singular than the classical solution. However those are the 
singularities of the field in the complex plane of t (or equivalently of z) which give 
rise to the factorial growth of the amplitudes. The appearance of the singularity at 
the imaginary axis follows from the simple fact that on this axis the classical field 
equation 

^0 = m 2 + A0 3 (31) 

corresponds to the free fall in the inverted A</> 4 potential which takes a finite time for 
a finite staring value 0(0). It looks at least extremely unnatural that quantum effects 
would slow down this fall to the extent that the time of the fall would be infinite. For 
any finite time, however, the singularity of the field there will produce the factorial 
growth of the amplitudes. 

As a simple final remark, it can be mentioned that though the present calculation 
is done for the case of unbroken symmetry it looks quite straightforward to apply the 
same technique to the case of the spontaneously broken symmetry. 

I am thankful to Lowell Brown for sending to me the preprint of his paper and 
to Arkady Vainshtein for stimulating discussions. I also acknowledge an extensive 
use of the Mathematical^ software in doing the calculations and preparing the text 
of the present paper. This work is supported in part by the DOE grant DE-AC02- 
83ER40105. 
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Fig.l The structure of the classical field <po{t) (eq.(§)) in the complex t 
plane. Heavy dots indicate the poles of 0o- The vertical line going between 
the poles is the axis of real r on which the operator (|13D is real. 
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